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ABSTRACT 

We explore the origin of mid-infrared (mid-IR) dust extinction in all 20 nearby (z < 0.05) bona-fi.de 
Compton-thick (Ah > 1.5 x 10 24 cm~ 2 ) AGN with hard energy (E > 10 keV) X-ray spectral measure- 
ments. We accurately measure the silicate absorption features at A ~ 9.7(xm in archival low-resolution 
(R ~ 57-127) Spitzer Infrared Spectrograph (IRS) spectroscopy, and show that only a minority 
(w 45%) of nearby Compton-thick AGN have strong Si-absorption features (S9.7 = ln(fmt/ fobs) <L 0.5) 
which would indicate significant dust attenuation. The majority (~ 60%) are star-formation dom- 
inated (AGN:SB<0.5) at mid-IR wavelengths and lack the spectral signatures of AGN activity at 
optical wavelengths, most likely because the AGN emission-lines are optically-extinguished. Those 
Compton-thick AGN hosted in low-inclination angle galaxies exhibit a narrow-range in Si-absorption 
(£9.7 ~ 0-0.3), which is consistent with that predicted by clumpy-torus models. However, on the basis 
of the IR spectra and additional lines of evidence, we conclude that the dominant contribution to the 
observed mid-IR dust extinction is dust located in the host galaxy (i.e., due to disturbed morphologies; 
dust-lanes; galaxy inclination angles) and not necessarily a compact obscuring torus surrounding the 
central engine. 

Subject headings: galaxies: active - Infrared: galaxies - X-rays: galaxies 



1. INTRODUCTION 

The unified m odel of active galactic nuclei (AGN; e.g., 
I Antonuccll 119931) is crucial to our understanding of the 
growth and evolution of massive black holes and their 
host galaxies. A key issue in AGN physics is the nature 
of the obscuring medium surrounding the central engine. 

There are three competing models for the specific 
structure, geometry and composition of the obscuring 
material surrounding the central supermassive black 
hole (SMBH): (1) a uniform (smoot h/continuous) , 
heavi ly -obscu ring, sub-parsec scale torus (|Pier fc Krolikl 
119921 I1993D : (2) an extended torus yi elding 
moderate obscuration (IGranato fc Danesd 119941; 
Efstathiou fc Rowan-Robinsonl 119951 : IGranato et al.l 
1997); and (3) clu mpy 'torus' of many indivi dual opti- 
cally thick clouds (iNenkova et alJl2002L 120081) . Each of 
these theoretical models make clear predictions for the 
observed torus properties. The graphite/silicate dust 
contained within the torus is predicted to extinguish 
AGN emission at UV/optical wavelengths, and should 
be co-spatial with cool neutral gas whi ch readily absorbs 
the X -ray emission from the AGN (e.g..[Mushotzkv et al.l 
1993). However, this dust/gas- rich torus is also predicted 
to (isotropically) re-emit in the mid-infrared (mid-IR; 
A ~ 5-50^m) with a mid-IR spectral energy distribution 
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characterized by a powerlaw-like AGN continuum, 
superposed with silicate absorption/emission features at 
9.7/iTO and 18^*m. The strength of these dust features 
are expected to be dependent on the specific geometry 
and optical depth of th e torus (e.g., iFritz et ahl 120061 : 
ISchartmann et aT]|2008D . To distinguish between these 
torus models, we require well-constrained results from 
sensitive and high-quality observations. 

In nearby optical and radio-selected AGN, a weak cor- 
relation between gas column density (Ah) and silicate 
(Si) absorption strength has been obse r ved fr o m Spitzer 
mid-I R spectroscopy (e.g., iShi et all 120061 : IWu et all 
2009). In turn, this suggests that selection of sources 
with strong Si-absorption is a good method to find the 
most heavily obscured AGN, i.e., Compton-thick sources 
with Nh > 1.5xl0 24 cm~ 2 (e.g.. lGeorgantopoulos et ahl 
120111 ). To first-order, these results would appear to agree 
with those predicted by simple uniform torus models. 
However, there is growing evidence that not all of the 
dust-extinction can be attributed to an obscuring cen- 
tral torus. 

U sing ground-based h igh-spatial resolution photome- 
try, [GandhTeFaD (|2009l ) show that Compton-thick AGN 
do not require significant corrections for dust extinc- 
tion to their AGN-produced mid-IR continuum, sug- 
gesting that Si-absorption features in Compton-thick 
AGN are not being produced in their nuclear regions 
(see also IHonig et al.l 120101: lAsmus et all 1201 ID . Fur- 
thermore, those AGN with strong Si-absorption fea- 
tures are often found to be hoste d in high l y in- 
clined and/or merging galaxies (e.g., iDeo et ahl 120071 
2009), providing first-order evidence that significant 
dust attenuation can occur within the AGN host-galaxy 
(e.g., lAlonso-Herrero et al.l 120111 and references there- 
in). Moreover, dust in the host-galaxy appears to extin- 
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guish the optical e mission-line signatures in » 25 -50% of 
nearby AGN (e.g.. iGoulding fc Alexanderl [20091) . and in 
a minority of case s, even the mid-IR AG N emission-lines 
(e.g., NGC 4945: IGoulding et al.l MHh . This apparent 
ambiguity among recent observations and theory raises 
two fundamental questions: (1) is Si-absorption a com- 
mon feature amongst Compton-thick AGN, as predicted 
by a unified AGN model? and (2) does Si-absorption 
predominantly arise from dust in a central torus or from 
within the AGN host-galaxy? 

To address these questions we explore the origins of 
apparent high optical-depth dust in the sample of all 20 
nearby (z < 0.05) Compton-thick AGN which are un- 
ambiguously identified with high-energy X-ray spectral 
observations. We use Spitzer-IRS observations to de- 
termine the average mid-IR spectral energy distribution 
of Compton-thick AGN, and establish whether silicate 
absorption is a signature common to the most heavily- 
obscured AGN. In §2 we describe the sample of Compton- 
thick AGN, outline the data reduction methods of the 
Spitzer-IRS observations and present the mid-IR spec- 
tra for the Compton-thick AGN sample. In §3 we com- 
pare the absorbing gas column and the apparent dust- 
extinction levels towards the central AGN and present 
the physical implications of these in light of a clumpy 
torus model. Finally, in §4 we review our results, find- 
ing that the mid-IR Si-absorption features observed in 
Compton-thick AGN are produced primarily due to dust 
in the AGN host-galaxy and not within a central torus. 



COMPTON-THICK AGN SAMPLE 
ANALYSES 



DATA 



Compton-thick AGN (7V H > 1.5 X 10 24 cm~ 2 ) are 
the most heavily obscured class of AGN; by their very 
nature, they are extremely difficult to detect and re- 
main hidden in most X-ray surveys (e.g.. iNorman et all 
2001 lAlexander et al.ll20"08l [201lt iComastri et all 120111: 



Feruriio et al.l 1201 It IGilli et all 12011b ILuo et al.l l20l"Tlf . 

Conclusive identifications of Compton-thick AGN are 
made through spectroscopic X-ray observations per- 
formed at E > 10 kcV where the relatively unabsorbed 
high-energy emission can be directly detected. The cur- 
rent sensitivities of E > 10 keV observatories (e.g., 
Beppo-SAX, Swift, Suzaku, INTEGRAL) are substan- 
tially limited by high backgrounds, relatively small ef- 
fective areas and low spatial resolution. To date, only 
20 bona-fide Compton-thick AGN have been unambigu- 
ously ident ified in the Universe a t E > 10 keV (for a 
revi ew see iDella Ceca et al.l 120081 and refe rences there- 
in; lAwaki et all 120091: iBraito et al.l [20091) . These 20 
Compton-thick AGN are all nearby (z < 0.05) systems 
hosted in spiral (Hubble-type SO or later) or merging 
galaxies with intrinsic X-ray luminosities Lx, intrinsic ~ 
(0.06-200) x lO^ergs" 1 (E ~ 2-10 keV), and hence 
span the wide range of AGN power observed in local 
Seyfert systems. 

2.1. Spitzer-IRS Data Reduction 

All 20 bona-fide Compton-thick AGN identified to date 
have archival low-resolution Spitzer-IRS spectroscopy. 
Specifically, these 20 AGN have been observed with the 
low-resolution modules (short-low [SL; 5.2-14.5/jm] and 
long-low [LL; 14.0-38.0/im]; R « 57-127) in either star- 
ing or mapping mode as part of multiple programs, and 



hence these data form a heterogeneous, but still complete 
sample. Our sample of nearby Compton-thick AGN and 
their basic properties are given in Table 1. The mid-IR 
data for many of our AGN sample have been analyzed 
using different reduct i on techniques in previous papers 
/e.g.. iDeoet^ 120071 IHao et all [20071: IWu et al.l l200l 
Mull anev et al.l 120101 |20TTT Sazonov et al. 2012, sub- 
mitted). However, in order to ensure self-consistency for 
the measurement of the spectral absorption/emission fea- 
tures which are integral to the analyses presented here, 
we have re-extracted each of the IRS observ ations us- 
ing our own custom IDL reduction routine (see | Gouldingl 
[20ia[Gmiiding et al.ll20TTb IMullanev et al.ll2011l) .~ 

Briefly, the two-dimensional Basic Calibrated Data 
(BCDs) images, produced by the S18.18.0 Spitzer Sci- 
ence Center (SSC) pipeline, were retrieved and rigorously 
cleaned of rogue 'hot' pixels using our customized ver- 
sion of IRSCLEAN. Next, individual rows were fit as a 
function of time to remove latent charge which exists on 
the detector after bright and/or long observations. The 
IRS Staring observations were averaged in the different 
nod positions, which were then used to perform alternate 
background subtractions of the source in each nod posi- 
tion and final spectra were extracted using the Spitzer- 
IRS Custom Extraction (spice) software provided by the 
SSC. Cleaned and processed IRS Mapping observations 
were input to the SSC package cubism to build the spec- 
tral data cubes and extract the 1-D spectroscopy. In- 
dividual IRS orders in the 1-D spectra were clipped of 
noise (see the Spitzer-IRS handbook for further informa- 
tion) and stitched together by fitting low-order polyno- 
mials to produce the final continuous spectra for each 
source. In the upper-panel of Fig. [T] we show the re- 
sulting low-resolution SpitzerARS spectroscopy for our 
complete sample of nearby Compton-thick AGN. 

2.2. Spitzer-IRS Spectral Decomposition & Silicate 
Absorption Measurement 

On the basis of uniform dust torus models 
produced from radiative tr ansfer theory (e.g., 
iDullemond &; van Bemmeil 120051 : Schartmann et al. 
2005; Fritz et al. 2006), the mid-IR spectral en- 
ergy distributions of intrinsically heavily-absorbed 
Type-2 AGN are expected to be dominated by an 
AGN-produced powerlaw and most, if not all, are 
expected to exhibit significant Si-absorption features 
at A ~ 9.7/iTO. Throughout this manuscript, we de- 
fine 59.7 as the depth of the Si-absorption feature at 

A — 9.7/OT7 where 59.7 = Zra(/g.74 n t r insic//9.7,obscrvcd). 

Following many previous studies, we use £9.7 as a 
good proxy for the apparent optic al depth and dust- 
extinction in our sample (e.g.. | Sppon et al.l 1 2007 : 
Shi et all 12001 iLevenson et al.l 12007 b iDeo et al.l 12007 . 
20091 : iGeorgantopoulos et al.l 120111 ) Though we note 



that the relation between S9.7 and true optical depth is 
most-likely non-line ar and somewhat model and orienta- 
tion dependent (see ILev enson et aTll2007t iNenkova et all 
2008; iSchartmann et al.ll2008h ~ 

Previous investigations have often employed simple ex- 
trapolation methods to estimate the depth of the silicate 
feature in mid-IR spectroscopy (e.g.. iSpoon et al.l [20071 : 
ILevenson et al.ll2007t IGeorgantopoulos et alll2011l) . This 
involves measurement of the observed continuum blue- 
ward and red-ward of the silicate feature, extrapolating 
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TABLE 1 
The Compton-thick AGN sample 
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Notes:- 'Conimon galaxy name: 2 J2000 positional co-ordinates from the NASA Extragalactic Database (NED); 3 Spectroscopic redshift; 4 Luminosity distance in Mega- 
parsccs assuming Hq = 71 kms -1 Mpc -1 and S"2 A = 0.70 and corrected for non-cosmological flows; 5_6 Galaxy morphology and major/minor axis ratio taken from the 
Third Reference Catalog of Bright Galaxies ( do VaupnuloarD ot ml^ ^Wj,) ; 7 Y _ denotes dust-lanes or peculiar morphology are evident in the optical image of the source; 
8 Optical classification in NED; 9 Logarithm of intrinsic (2-10 keV) X-ray luminosity in units of orgs" 1 corrected for Compton-thick absorption for those sources with 
JV H < 10 25 cm -2 ; l°Ncutral hydrogen column density derived from hard X-ray spectroscopy in units of 10 24 cm~ 2 ; U Rcfcrcncc for X-ray data (sec below); 12 Typc of 
.S'pifzer-IRS observation (M: Mapping; S; Staring); i3 Spitzer observation record indicator; l 4 Datc of Spitzer-IKS observation; l^Dcpth of Si-absorption at A ~ 9.7ttm derived 
from spectral fits to Spitzer-IRS data using the formalism of Mullaney ct al. (2011); l 6 Depth of Si-absorption at A ~ 0.7fim derived from spectral fits to Spitzcr-IRS data 
using PAHFIT; 17 AGN-starburst ratio derived from fitting IRS spectrum with DecompIR (AGN-dominated system = 1). 

References:- (1) ilwmmra ot al. ( 3 001, ); (2) . Collirmo fa BiDn dt 1 3000 .1 ; (3) . ■■, 1 . . 1 I ot al. <-UWtJ, i: (4) . Snnonov ot ol„ ( 3007 , ); (5) .Boolimonn ot al . ( 3006 , ); (6) iBnooani ot al. 
(2007); (7) Loto-sco.!-, ot al. (2006); (8) Ucda ot al. (2007); (9) Cappi ot al. (1900); ( 1 0) Mark-.vardt ot al. (2005); (11) Baooani ot al. (1900); (12) Av.-aki ot al. (2000); (13) 
I .,,, nn , I , ot al. WQ8J. 1; (14) . VimnU ti Cmim lri WW ); (15) . Mainline ot oh . (100 & 1; (16) , Polio Coaa ot a t. (300 a) ; (17) . Hallo at nl .. (3884.1 ; (18) . Benito ot al . (3084, 1; (19) 

C ol 0,1 ., (3888.) ; (20) . Fulrannwa, ot 0,1 . ( 3 001 .1 ; (21) . Mot* ot ol .. t+OOft) ; (221 . Iwnonwa nt al. ( 1007 , 1; (231 . Comootri ot nl . (3007 . 1; (241 . Viimoti ot ol ., (1000 . 1; (25) 

Turner et al. (2000); (26) Malaguti et al. (1998); (27) Braite ot al. (2009) 



a powerlaw between the continuum points and compar- 
ing this to the observed flux at A ~ 9.7/im. In prin- 
ciple, this extrapolation method is sufficient for AGN- 
dominated spectra. However, as shown in Fig. [TJ the 
majority of the mid-IR spectra for Compton-thick AGN 
contain strong polycyclic-aromatic hydrocarbon (PAH) 
features (at 6.3, 7.7, 11.3 and 12.3/im) indicative of sub- 
stantial circumnuclear star-formation activity. These 
PAH features provide significant contributions to the 
mid-IR spectra blue-ward of the Si-absorption feature, 
having the effect of artificially increasing the continuum 
strength at A ~ 6-8fj,m, and hence, over-estimating the 
intrinsic AGN continuum flux at A ^ 9.7/im. In a pure 
synthetic star-forming template with no evidence for ex- 
tinction, a substantial Si-absorption feature (S9.7 ~ 0.7) 
would be inferred using this basic method. Therefore, 
we choose to use a spectral decomposition method to re- 
move the host-galaxy emission, and measure the depth 
of the Si-absorption feature only within the AGN com- 
ponent. Furthermore, to constrain any systematic uncer- 
tainties derived from using any one spectral decomposi- 
tion method we use two independent spectral decompo- 
sition procedures. 

First, we use the IDL software package DecompIR (Mul- 
laney et al. 2011), which performs a Chi-squared mini- 
mization of the observed spectrum to a combination of 
starburst templates (i.e., a host-galaxy component) and 
an absorbed power-law (i.e., an AGN comp onent) . Here , 
we employ a simple screen extinction curve (|Drainell2003| ) 
to account for absorption of the AGN component; thus, 
we do not attempt to constrain the physical region where 
dust extinction is occurring. Following Goulding et al. 



(2011), we allow for a range of star-forming templates, 
including that of the archetypa l nearby starburst, M82; a 
combined Brandl et al. (2006]) starburst template^ and 
a range of compact nucl ear and extra-nuclear theoret i- 
cal starburst templates (jSiebenmorgen fc Kriigell [20071, . 
From the fitted AGN component, we predict the intrinsic 
power of the AGN at A ~ 9.7/im, and compare this to 
the observed luminosity to establish an estimate for S9.7 
(column 15 of Table 1). Second, we employ the widely 
used spectral-fitting package, PAHFIT (jSmith et al.1120071) 
to fit a range of heated dust continua combined with 
PAH and emission line features to derive 6*9.7 (colum n 
16 of Table 1). Following. iGeorgantopoulos et ali (|2011f) . 
we change the default dust continuum temperatures in 
PAHFIT and additionally allow hotter dust temperatures 
(400-1400K) to be fit to the spectra, these are expected 
to arise from an accretion heated AGN torus. In gen- 
eral, we find very good agreement for S9.7 derived using 
DecompIR and PAHFIT (see columns 15 and 16 of Table 
1). In U21 we compare these two measures of the Si- 
absorption feature to gas column density in Compton- 
thick AGN. 

2.3. The Average mid-IR Spectra of Compton-thick 

AGN 

In Figure[T]we show that Compton-thick AGN at z ~ 
are characterized by a variety of spectral shapes with 
AGN-starburst ratios of - 0.1-1. Of the 20 AGN in 
our sample, only four appear to be dominated by a fea- 

6 We do not include NGC 660, NGC 1365, NGC 3628 and 
NGC 4945 as part of the Brandl et al. (2006) template as these 
sources arc known to harbor AGN. 
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Fig. 1. — Upper Panel: Low resolution (i?, ~ 57-127) Spitzer- 
IRS spectroscopy of all 20 bona-fidc Compton-thick AGN in the 
nearby Universe with high-energy (E > 10 keV) spectroscopic con- 
straints. The gray shaded region represents the wavelength bound- 
aries for the dominant Si-absorption feature (A ~ 9.7fim) in the 
short-low (SL) Spitzer-TRS pass-band. The mid-IR spectral energy 
distributions of Compton-thick AGN are clearly very diverse with 
spectra dominated by a combination of an AGN-produced power- 
law; strong starburst-produced polycyclic aromatic hydrocarbon 
features at A ~ 6.3, 7.7, 11.3 and 12.8/jm; and in a minority 
(9/20) of cases, strong (S9.7 > 0.3) Si-absorption features. Ad- 
ditionally, we show the average mid-IR spectrum of the sample of 
Compton-thick AGN (thick black line). Lower panel: Spectral 
decomposition of the average Compton-thick AGN spectrum (blue 
solid curve) The best-fit extinction-convolved AGN power-law and 
starburst template are shown with dotted and dashed curves, re- 
spectively. The total best-fit spectrum (i.e., power-law + starburst 
+ emission-lines) is shown with a solid red curve. The best-fitting 
constraints for various AGN classifications are given in Table 2. 



tureless power-law continuum, while eleven show strong 
PAH features indicative of circumnuclear starburst activ- 
ity which, for seven (w 35%) of the sources, even dom- 
inates the bolometric output of the galaxies at mid-IR 
wavelengths (AGN:SB < 0.5; see Column 17 of Table 
1). By contrast, Sazonov et al. (2012, submitted) find 



that only < 20% of all hard X-ray selected Compton-thin 
(AT H ~ (0.05-100) x 10 22 cm" 2 ) Type-2 AGN detected by 
INTEGRAL are dominated by star-formation at mid-IR 
wavelengths. It is therefore possible that enhanced star- 
formation activity may be a characteristic specific to the 
most heavily absorbed (A H > 10 24 cm~ 2 ) AGN, and not 
merely a hard X-ray selection effect, at least in the nearby 
Universe. 

Through spectral stacking, we find that Compton-thick 
AGN, on average, show strong star-formation activity. 
We stacked all 20 mid-IR spectra to produce a mean mid- 
IR spectral energy distribution for Compton-thick AGN. 
We chose to normalize the rest-frame mid-IR spectra for 
each of the Compton-thick AGN at A ~ 6.5-7 /im (i.e., 
a wavelength range with a relatively featureless contin- 
uum). In the lower panel of Fig. [2 we show the spectral 
decomposition of the av erage spectrum mod eled using 
DecompIR. Our modified iBrandl et al.l J2006) template 
provides the best x 2_nt 'starburst' to the mean Compton- 
thick AGN spectrum, and contributes « 40% of the over- 
all luminosity to Compton-thick AGN at mid-IR wave- 
lengths. 

As we show in Table 1 and Fig. [TJ the detection of Si- 
absorption features in nearby Compton-thick AGN is far 
from ubiquitous. The Compton-thick AGN cover a wide 
range in Si-absorption depth, 6*9.7 ~ 0-1.4. While the 
mid-IR spectra for eight of the Compton-thick AGN do 
appear to be moderately extinguished, we also find that 
the spectral decomposition (using DecompIR) of NGC 424 
shows strong evidence for Si-emission (S9.7 < 0.0); based 
on unified models, Si-emission features are expected only 
in unobscured Type-1 AGN, irrespective of specific torus 
geometry /composition. This result is also partially con- 
firmed by PAHFIT which produces a relatively poor fit 
(Xn-i ~ 2 6 - 1 ) for NGC 424 and b y design does not al- 
low apparent optical depths of Sg 7 < 0.0. Explanations 
for Si-emission features in Type-2 A GN (e.g., NGC 2110 
Mason et al. 2009; SST1721+6012. INikutta et^[2f)M ) 
are that the emission is either (i) arising from the inner- 
most region of the AG N narrow-line regio n, above the 
scale-height of the torus ()Mason et al.ll200"9l ); or (ii) from 
an unobscured inner region of a low-mass clumpy torus 
(e.g., based on the models of Schartmann et al. 2008). 

From the AGN powerlaw component of the stacked 
Compton-thick AGN spectrum, we directly measure the 
depth of the Si-absorption feature at A ~ 9.7/im. We 
find that, on average, the re-radiated AGN emission of a 
Compton-thick source experiences an apparent attenua- 
tion of only S9.7 ~ 0.36±0.04, equivalent to an optical ex- 
tinction of Ay w 3.9-4.9 magsQ Optical depths such as 
these are readily observed in nearby galactic star-forming 
regions (e.g., Draine 2003). In § 13.21 we further explore 
whether this measured dust extinction in Compton-thick 
AGN can be attributed only to that arising in the host- 
galaxy and not necessarily from an obscuring torus. 

3. GAS ABSORPTION AND DUST EXTINCTION 
IN COMPTON-THICK AGN 

Unified AGN models postulate that obscuring gas 
and dust are both located in the circumnuclear torus, 

7 Throughout w e assume a standard optical-IR extinction curve 
HDraine &; Lil2007h in the direction of the Milky- way galactic center 
such that Ag.r,j. m /A v ~ 0.075. 
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Fig. 2. — Si-absorption depth at 9.7fim (Sg.7; a good proxy for apparent dust extinction) observed in Spitzer-IRS mid-IR spectroscopy 
versus hard X-ray derived gas column density (A^h) in units of cm -2 for all bona-fide Compton-thick AGN in the nearby Universe. Error 
bars represent the range in Sg.7 measured using DecompIR and PAHFIT. Constant gas— dust ratios are shown (dotted lines) for ratios of 
A^h,x/A ? h,ir — 100 and 1000 assuming Nn/Ay = 1.93 X 10 21 cm" 2 and Ag. 7 /A v ~ 0.075 (e.g,. IDraine fc Lj [2007). Inclination angles are 
shown with color-gradient and are derived from major-minor axis ratios (see Column 6 of Table 1), face-on systems (b/a ~ 1) are shown 
in blue, edge-on (b/a ~ 0) in yellow, colors in between are linearly scaled to intermediate inclination angles (b/a ~ 0-1), and merging 
systems (established from optical imaging and NED classifications) arc shown in red. Optical spectral classifications (AGN; LINER; 
HH/star- forming) are shown with filled circles, triangles and stars, respectively. Additionally, sources with dust-lanes or disturbed optical 
morphologies are highlighted with open squares (see Column 7 of Table 1). For comparison, we show the linear-fit presented in Shi et al. 
(2006) for a sample of Seyfert galaxies with measured ATjj. 



so that large gas columns are associated with heavy 
dust extinction and, for a uniform torus, deep Si- 
absorption features. With this in mind, a number of 
recent studies use the existence of high apparent optical 
depth (5*9.7 > 1) to identify C ompton-thick AGN (e.g. 



Georgantopoulos et all 120111 : INardini fc Risalitil 120111 
Dasyra et alJ 12011). However, as shown in Fig. 1, the 
large range in S9.7 seen across our complete Compton- 
thick AGN sample suggests that the dust producing the 
Si-absorption features may not be co-spatial with the X- 
ray absorbing gas. 



3.1. Compton-thick AGN with little apparent dust 
extinction 

To establish if a physical link exists between the gas 
and dust in Compton-thick AGN, as expected by a uni- 
fied model, in this section, we compare the depth of the 
Si-absorption feature and the gas column density. In 
Fig. [21 we plot the measured X-ray gas column density 
versus mid-IR Si-absorption strength (59.7). For each 
AGN, we use the average value for Sg.7 derived inde- 
pendently from DecompIR and PAHFIT in the previous 
section. Error bars represent the spread in the values of 



6 



A.D. GOULDING ET AL. 



5*9.7 for a particular source. We find little or no observed 
correlation between gas absorption and dust extinction 
in Compton-thick AGN (Spearman's rank, p ~ —0.30, 
-fnuii ~ 0.21). At least eight of the sources exhibit 
little or no evidence for dust extinction (S9.7 < 0.1), 
and moreover, NGC 424 has silicate emission features 
(5*9.7 < 0.0), in general, such emission features require 
a direct line-of-sight to the inner AGN torus. By con- 
trast, 9 of the 20 (~ 45 ± 18%) Compton-thick AGN 
in our sample appear moderately (6*9.7 ~ 0.5-1.0) to 
heavily (5*9.7 > 1-0) obscured at mid-IR wavelengths, as 
would be expected by unified AGN models that invoke 
a smooth uniform dust torus. However, the fraction of 
Compton-thick AGN with high apparent optical depth 
(69.7 ^ 0-5) is similar to that measured for lower column 
density {N H ~ (0.05-50) x 10 22 cur 2 ) AGN detected 
by Swift-BA T, 30 ± 10% (derived from the spectroscopy 
presented in iMullanev et al.ll2010[ ). suggesting that large 
Nh does not necessarily predicate large dust extinction. 
Indeed, at least fo ur (sa 10%) Compton- thin Swift-BAT 
AGN presented in IMullanev et all ((20101 ) have S9.7 > 1. 

In Fig. [5] we additionally show curves of constant gas 
column density predicted using 69.7 assumi ng a standard 
optic al-IR Milky-way extinction curve ([Draine fc Lil 
12007ft and a gas -dust ratio of Nn/A v = 1.9 3 x 
10 21 cm" 2 mag" 1 ([Fitzpatrickl [19851: iDraind [20031 for 
R v = A V /E{B-V) = 3). We show that Compton-thick 
AGN are required to have gas-dust ratios which are a fac- 
tor ps 50-10,000 greater than observed galactically, which 
may suggest there is no single g as-dust ratio in these 
sources (e.g. lMaiolino et aT1l2001l ). These large gas-dust 
ratios may be explained by additional absorption of X- 
ray emission by large quantities of neutral gas residing 
within the dust-sublimation radius (0.1 < r < 10 pc) of 
the obscuring torus, so-called 'partial covering'. 

Recent X-ray monitoring campaigns of variable absorp- 
tion state AGN (e.g., Mrk766; NGC 1365; NGC 4151) 
have observed short-term occultations of the X-ray 
source by high- velocity (~ 3000 km s" 1 ) circumnuclear 
gas cl ou ds residing at ~ 0.1 pc (e.g., IRisaliti et al.l 
I2009allb1 : iWang et all 120101 ). The motion and distance 
from the central source of these gas clouds is consis- 
tent with that expected for X-ray abso rption arising 
from the AG N broad- line region (BLR; e.g. . IRisaliti et all 
I2009U 1201 If) . Within these relatively small radii, SiO 
dust grains are preferentially destroyed due to the high 
temperatures (T > 1000K). Hence, the gas-dust ratio 
of the BLR may be far larger than that of the dusty- 
torus. However, partial-covering may not explain the 
high gas-dust ratio observed in all of the Compton-thick 
AGN included here. R ecent, spatially-resolved results of 
iMarinucci et al.l (|2012f l suggest that the majority of the 
Compton-thick absorption in NGC 4945 may be at far 
larger radii (~ 100 pc) than the BLR. These variations in 
BLR and torus absorption s, as well as inclination effects 
(e.g., iNenkova et aT] 12008) may add substantial scatter 
to an expected gas-dust ratio for nearby AGN. Taken 
together, these results give rise to several interesting sce- 
narios for: (1) the spatial extent of the intervening gas 
and dust; (2) the density gradient of the materials across 
the torus; and (3) possible gas-dust separation bound- 
aries within a torus. Given our small sample and (by 
design) our limited range in N-r, these issues are beyond 



TABLE 2 

Stacked properties of Compton-thick AGN 



Sub-sample 


i'9.7 


AGN:SB 


(1) 


(2) 


(3) 


All sources 


0.36 


0.62 


Optical Seyfert* 


0.14 


0.83 


Optical LINER/HII* 


0.55 


0.43 


Face-on/intermediate Inclination 


0.19 


0.75 


Edge-on Inclination or Merger 


0.96 


0.44 



Notes:- 1 Host galaxy properties of sub-sample of main Compton- 
thick AGN sample; 2 Average depth of Si-absorption feature (S9.7) 
within stacked mid-IR spectra of sub-sample, measured using 
DecompIR; 3 AGN-starburst ratio derived from decomposition of 
stacked mid-IR spectra using DecompIR (AGN-dominated = 1). 
* Excluding QSOs. 

the scope of this investigation. We may still robustly 
conclude that all Compton-thick AGN (i.e., large JVh 
sources) are not necessarily obscured by cool dust with 
high apparent optical depth. However, the reverse of this 
statement may still also be true, in that all systems ob- 
served to have 69.7 3> 1 may all contain Compton-thick 
AGN, but are currently missed due to the lack of suffi- 
ciently sensitive instrumentation; although these extreme 
69.7 3> 1 sources are only a limited subset (< 20%) of 
the Compton-thick AGN population. 

3.2. Evidence for dust extinction arising in the host 
galaxy of Compton-thick AGN 

The apparent lack of a correlation between gas and 
dust attenuation in Compton-thick AGN presented in 
the previous section sugge sts that the gas and dust may 
not be co-spatial. Indeed, ISturm et al.l (|2005D suggested 
an extranuclear origin for the Si-absorption features ob- 
served in both low and high-luminosity AGN. This is 
supported by Gandhi et al. (2009) showing in high spa- 
tial resolution ground-based photometry of nuclear re- 
gions in nearby AGN, that the AGN-produced mid-IR 
continuum appears almost isotropic and unaffected by 
Si-absorption. Furthermore, the inclination angle of the 
host galaxy may have a significant obscuring effect on 
the AGN emission (e.g.. iMalkan et al.lll998t iMatt et aT] 
I2000t lAlonso-Herrero et al.l 120111: lLagos et all 120111) . In 
particular, those AGN presenting strong Si-absorption 
features are often found t o be hosted in highly-in clined 
and/or merging galaxies ([Deo et alJ 120071 12009). For 
some host-galaxy types and inclinations, the dust extinc- 
tion may be so over-whelming that it has the ability to 
completely obscure an AGN (e.g., at optical and even 
mid-IR wavelengths; IGoulding fc Alexanderl [20091) . In- 
deed, the mid-IR AGN emission in the late-type edge-on 
AGN, NGC 4945 is almost entirely obscured based on 
Spitzer-IRS spectroscopy. Hence, in this section, we test 
for a link between the observed mid-IR extinction and the 
properties of the host galaxy in Compton-thick AGN. 

In Fig. [2j we highlight the optical classifications and 
host-galaxy inclinations for the Compton-thick AGN 
sample^ As noted previously, each of the 20 Compton- 
thick AGN are hosted in late-type (SO-Scd) or irreg- 
ular/merging galaxies. From stacked mid-IR spectra, 

8 The host-galaxy inclination angles are derived from the major- 
minor axi s ratio as defined in the Thir d Reference Catalog of Bright 
Galaxies (de Vaucoulcurs et al. 1991). 
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we find that the galaxies classified optically as Seyferts 
(i.e., using standard optical emission- line diagnostic di- 
agrams) have lower average apparent dust extinction 
(5*9.7 ~ 0.14) than the Compton-thick AGN hosted in 
the optically-classified Low-Ionization Nuclear Emission 
Line Regions (LINERs; i.e., sources with optical spec- 
tra characteristic of LINERs) and HII galaxies (6*9.7 ~ 
0.55). This dichotomy between apparent dust extinc- 
tion and galaxy classification provides further evidence 
that the optical class may be influenced by high-levels 
of intervening dust in the host galaxy. Additionally, we 
show that Compton-thick AGN hosted in face-on (low- 
inclination) galaxies have mid-IR spectra exhibiting lit- 
tle or no apparent dust extinction (S9.7 ~ 0.19). By 
contrast, Compton-thick AGN hosted in edge-on and/or 
merging systems have significantly greater dust extinc- 
tion (5*9.7 ~ 0.96; see Table 2) clearly suggesting that the 
Si-absorption features observed at mid-IR wavelengths 
have an extra-nuclear origin. 




Wavelength (jjm) 



Fig. 3. — Schematic diagram illustrating the effect of host-galaxy 
inclination (and morphology disturbance) on the perceived dust 
extinction at mid-IR wavelengths. Based on the effects presented 
here, the strength of the Si-absorption feature at A ~ 9.7fj,m qual- 
itatively appears to increase with increasing galaxy inclination, 
i.e., consistent with the view that the obscuring dust which dom- 
inates the mid-IR extinction features arises only from within the 
host galaxy. Image credits, from top: NGC 5194'. NASA and 
the Hubble Heritage Team (STSci/AURA); NGC 4939: Adam 
Block/NOAO/AURA/NSF; NGC 4945: ESO; Arp 299: NASA 
and the Hubble Heritage Team (STSci/AURA) 

There are also three optically-classified Seyferts 
(NGC 3281; NGC 4939; Circinus) within our Compton- 
thick sample which have larger Si-absorption features 



(59.7 ~ 0.8-1.5) than would necessarily be expected 
given their moderate inclination angles (b/a ~ 0.5). 
Visual inspection of optical images used to derive the 
morphological classifications shows that these three 
Compton-thick AGN have disturbed nuclear morpholo- 
gies and/or dust-lanes aligned along the galaxy nucleus. 
These additional extra-nuclear features may provide fur- 
ther regions of cool-dust, and hence, larger Si-absorption 
strengths. Taken together, these results suggest that the 
dominant source of dust-extinction observed at mid-IR 
wavelengths is not necessarily linked with an obscuring 
torus, but instead is arising from extended kilo-parsec- 
scale regions which lie far outside the expected dust- 
sublimation region of an accreting SMBH. Indeed, this 
galactic scale origin for the cool dust, which is required 
to produce deep Si-absorption features, is consistent with 
the wide range of 69.7 observed in near by starburst- 
dominated Ultra-luminous IR galaxies (e.g . , iSpoon et al.1 
2007; lDesai etaLl [20071: iFarrah et al.ll2008l) . 

Despite our relatively small sample of sources, the 20 
Compton-thick AGN considered here cover a wide range 
of apparent dust extinction values, 5g.7 ~ 0.0-1.7. How- 
ever, in Fig. [2] we note a clear dearth of sources with 
5 9 . 7 > 0.3 and A H > 5 x 10 24 cm~ 2 . AGN in this regime 
are likely to be heavily obscured at both mid-IR and op- 
tical wavelengths due, in part, to large quantities of inter- 
vening dust within the host galaxy, as well as being heav- 
ily obscured at X-ray energies due to large gas columns 
within the torus. Such heavily-obscured AGN are still 
not robustly identified in the nearby Universe. Obser- 
vational expectations are that the number of Compton- 
thick AGN with A H ~ (1-5) x 10 24 c m" 2 and JV H >5x 
10 24 cm~ 2 are roughly co nstant (e.g.. lRisaliti et alll!999t 
iSalvati fc Maiolino|[2000]) . Based on the 5 9 . 7 -A H distri- 
bution of our Compton-thick AGN sample, we should 
expect w 6-8 sources in this 'heavily-obscured' parame- 
ter space; however, we find only one source, suggesting 
> 80% still remain unidentified in the nearby Universe. 
From the sources currently identified with 5g.7 > 0.5, we 
expect the majority of the 'heavily-obscured' systems to 
be optically-unidentified due to large quantities of host- 
galaxy dust as well as being starburst dominated in the 
mid-IR. Hence, they will still remain unidentified in high- 
redshift multi-wavelength blank-field surveys. This will 
potentially cause a bias in our current perceived view of 
the heavily-obscured AGN population and their place in 
galaxy evolution models. 

As we have shown, the Si-absorption features in many 
edge-on Compton-thick AGN appear to be unassociated 
with the dusty torus surrounding the central engine; we 
therefore cannot use these sources to constrain theoret- 
ical models for the AGN infrared emission. However, 
based on our results those Compton-thick AGN hosted in 
face-on galaxies appear relatively unobscured in the mid- 
IR with 5 9 . 7 ~ 0-0.3. Such AGN are unlikely to have 
significant host-galaxy dust extinction along the line-of- 
sight; hence, any observed Si-absorption will be intrinsic 
to the torus. For a sensible (physical) range of AGN 
torus parameters (see Table 1 of Fritz et al. 2006; Table 
1 of Schartmann et al. 2008), clumpy torus models (e.g., 
Nenkova et al. 2002; Schartmann et al. 2008), where-by 
the dust is randomly distributed in discrete dust clouds, 
predict similarly weak Si-absorption (5g.7 < 0.3) for 
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optically-thick edge-on (i = 90°) tori. By contrast, uni- 
form torus models (e.g., Pier & Krolik 1992; Fritz et al. 
2006) predict large silicate features (S9.7 3> 1) for similar 
input parameters. Indeed, for a Compton-thick system, 
uniform torus models reproduce Sg.7 ~ only when in- 
voking a steep dust-density gradient across an extremely 
compact (Rmax/Rmin "C 30) and somewhat unphysical 
torus. Therefore, under the paradigm of a clumpy torus, 
we suggest that the additional Si- absorption measured in 
those Compton-thick AGN which are hosted in edge-on 
host-galaxies becomes entirely consistent with extinction 
arising from dust which is distributed at very large radii 
(3> parsecs) within the main host galaxy. In light of these 
findings, in Fig. 3 we present a schematic cartoon of the 
effect of the host-galaxy inclination and morphology on 
the observed mid-IR Si-absorption feature. 

4. SUMMARY & CONCLUSIONS 

We have investigated mid-infrared (mid-IR) dust ex- 
tinction in all hard X-ray (E > 10 keV) detected, bona- 
fide Compton-thick (JVh > 1.5 x 10 24 cm" 2 ) AGN in 
the nearby Universe to constrain the dominant source 
of attenuating dust in heavily obscured AGN and build 
a more complete understanding of Compton-thick AGN 
at mid-IR wavelengths. We measured Si-absorption fea- 
tures at 9.7/im in archival low-resolution (R ~ 57-127) 
mid-IR Spitzer-IRS spectroscopy to infer the attenuating 
dust column in these highly gas absorbed systems and de- 
termine the average mid-IR properties of Compton-thick 
AGN. We found that the mid-IR SEDs of Compton-thick 
AGN are diverse, and only in a minority of sources (9/20) 
do we find strong Si-absorption features (6*9.7 ^ 0.5). In 
turn, we suggest that sample selection based solely on 
the mid-IR absorption features (for S9.7 > 1) may miss 
> 80% of Compton-thick AGN. From analysis of the 
average mid-IR spectrum, we find that Compton-thick 
AGN are characterized by roughly equal (40:60 ratio) 
quantities of starburst and AGN components which is 
only mildly dust obscured (S9.7 ~ 0.36 ±0.04; A w « 3.9- 
4.9), i.e., average dust-extinction levels which are con- 
sistent with those which are observed in host galaxy 
star-forming regions. Furthermore, the high-levels of 
star-formation observed in these Compton-thick AGN 
appears to be a feature specific to the most heavily- 
obscured AGN in the nearby Universe. If such a trend 
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continues towards higher redshifts, then infrared photo- 
metric and X-ray surveys which focus on the detection 
of powerlaw components to infer AGN actiyity w i thin a 
source (e.g., iDaddTetanHQQl IDonlev et al]l2007l l2008h 
will be partially biased against the detection of a signifi- 
cant proportion of the Compton-thick AGN population, 
and hence, may still miss the majority of the most heav- 
ily obscured AGN. 

We compared Nh and 6*9.7 m Compton-thick AGN, 
and find no significant correlation between these two 
measures of apparent obscuration. We find that the 
most heavily attenuated sources (S9.7 > 0.5) appear to 
be hosted in galaxies with visible dust-lanes, disturbed 
morphologies and/or galaxies which are highly-inclined 
along the line-of-sight. By contrast, we show that sources 
hosted in face-on galaxies present only weak dust extinc- 
tion features (S9.7 ~ 0-0.3). A similarly narrow-range 
in Si-absorption is predicted by those theoretical torus 
models which invoke clumpy distributions for the ob- 
scuring material. We suggest that the deeper silicate 
features (S9.7 > 0.5), which are observed in only a mi- 
nority of Compton-thick AGN (w 40%), arise from in- 
tervening dust at much larger scales than predicted for 
a torus (i.e., from within the host galaxy). When com- 
bined with previous investigations, this provides further 
evidence that the obscuring cool dust, which dominates 
the extinction seen at mid-IR (and optical) wavelengths, 
may not necessarily be co-spatial with a gas-rich central 
torus. 
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